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We introduce a D-dimensional Hamiltonian formalism for the study of Polyakov loop models of 
finite temperature gauge theories in D+l dimensions. Polyakov loop string tensions are obtained 
from energy eigenstates of the Hamiltonian. For D=l, the gauge theory reduces to quantum 
mechanics on the gauge group; for D>1, the Hamiltonian includes hopping terms that link sites on 
the transverse lattice. The deconfmed phase is associated with a ground state which breaks Z(N) 
symmetry, and Svetitsky-Yaffe critical universality emerges naturally for D>1. A minimal model 
is proposed which naturally reproduces approximate Casimir scaling for a range of couplings. 
Different classes of potentials lead to different pictures of how confinement is realized. Such 
potential energy terms also modify string tension scaling laws, as we demonstrate using two 
potentials: one representing the perturbative thermal contributions from gluons, and the other 
arising from magnetic monopoles in certain confining supersymmetric theories. 
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1. Introduction 

There is a substantial body of evidence from lattice simulations [0? 0> Ml B & 01 tnat str i n S 
tensions in confining gauge theories obey scaling laws. One strong possibility is Casimir scaling, 
where the ratio of string tensions obey G R > /or = Czri/Cir with C2R the quadratic Casimir for the 
representation R. Another possible scaling law, suggested by string theory, is sine-law scaling, 
in which the lowest string tension in each Af-ality sector obeys O^/Oy = sm(lc7t/N) /sin (ft/Af) 
[§, ^]. We will explore string tension scaling using a class of Hamiltonian effective models for 
Polyakov loops obtained from underlying D + 1 -dimensional gauge theories using the same type of 
arguments which lead to the Euclidean effective action for Polyakov loops [[K]]. These models are 
defined on a D-dimensional space with one continuous dimension and D — 1 discrete directions. 
The general form of the Hamiltonian may be written as 

H = Z i KC V + Vj] - % £ [XFjXh + XpjXFk] (1-1) 

J (Jk) 

The first term is the kinetic term for each site j, given by the Casimir operator. Higher-order 
terms may also be present. The second term is the potential on each site, Vj, a Z(N) -invariant 
class function of the group. The third term is a hopping term between nearest-neighbor sites, 
and generalizes to include other representations and longer hops. H acts on wave functionals 
which are class functions *P [P] of Polyakov loops. The Hilbert space is spanned by products 
of characters, with Haar measure providing the natual inner product. The integration over each 
Polyakov loop need only be taken over the maximal Abelian subalgebra, so that each Polyakov 
loop Pjm&y be represented as Pj = diag [e' dl ..e' 0N ~\ . The ground state energy density Eq in the 
Hamiltonian formalism is related to the free energy density / and the pressure p via f5f = —pp = 
E . 

A central question we address is naturalness: what class of potentials lead to Casimir scaling 
or similar behavior? Many potentials have regions of parameter space which tend to localize wave 
functions around N degenerate minima related by Z(N) symmetry. In such cases, Z(N) symmetry 
can be maintained in the confined phase by tunneling. If we have an approximate wave function 
<J>o (P) which is localized near one of the degenerate minima, we can construct a set of N wave 
functions as 

1 N-l r 



InijkfNQ ( e ~2itijk/Np 



(1.2) 



It is easy to show that the string tension splitting in this case is determined by tunneling, giving a 
scaling law of the form Ck/o\ = sin 2 (kn/N) / sin 2 (n/N). Another, complementary limit localizes 
the wave function around the unique set of eigenvalues left invariant by Z(N) symmetry. For SU (3), 
this gives the matrix Pq = diag [l,£ 27r '/ 3 ,e 4m / 3 ] . For any representation R with non-zero A^-ality, 
Xr (Pq) = 0, so a configuration of Polyakov loops of the form g(x)Pog~ l (x) can be said to confine. 
These limiting behaviors represent the extremes of the potentials we discuss below. 

2. Mean Field Theory 

In a Hamiltonian formalism, mean field theory is equivalent to the Hartree approximation. We 
assume that the ground wave function is a product of uncorrected site wave functions: [P] = 
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p=1/T 




Figure 1: Lattice geometry for Hamiltonian formalism. 

n Yj (Pj)- Minimizing the energy and requiring self-consistency gives the Hamiltonian Hj for the 
site j: 

Jrrd 

Hj = KC 2j + Vj — y~ [Xfj (Xf) + X*fj (Xf)] (2- D 

In the confined phase, Z(N) symmetry requires (xf) = 0, and we are left with a Hamiltonian for 
Z(A^)-invariant quantum mechniacs on the group manifold: 

[KC 2j + Vj] Win (Pj) = E n Yjn (Pj) (2.2) 

As a first approximation, we associate the eigenvalues with string tensions via E n — Eo = a n /T. 

The Hamiltonian formalism is capable of reproducing the behavior of the deconfinement 
phase transition revealed by familiar models based on an effective action |[T(]]. In the case of 
SU(2), assume for simplicity a potential Vj = —h A % A j and a trial wave function of the form 
Yoj = Vl — A 2 — B 2 + A%Aj + B%fj where A and B can be taken to be real. The order parame- 
ter for the deconfinement transition is the expected value of the Polyakov loop in the fundamental 
representation = m (A,B) = 2B (a + y/\ — A 2 ~B 2 ^ ■ The ground state energy density is given 
by the minimum of 

E var = k (C A A 2 +C F B 2 ) — (D— l)J F m 2 (A,B) - h A (A 2 + 2A + B 2 ) (2.3) 

The low-temperature, confined phase occurs when the term in E var proportional to K dominates, 
and £ vor is minimized when B = 0. If h A is non-zero, A will be non-zero, and (xa) W1 U have a 
non-zero expectation value at all temperatures. As the temperature grows, the hopping parameter 
Jf grows as well. When the coefficient of B 2 in E var becomes negative, there is a second-order 
phase transition to the deconfined phase where B ^ and Z(2) symmetry is spontaneously broken. 
Note that the h A potential term causes the critical temperature to change, but not the order of the 
transition. 

In the case of SU (3), we take a trial wave function of the form Yoj = \f\—A 2 — 2B*B +A% A j + 
B*XFj+B%pj and the order parameter is given by (XFj) = m (A,B) = 2B^A + y/\ —A 2 — 2B*B^j + 
(B*) 2 . The variational estimate of the ground state energy density is given by 

E var = k(C a A 2 + C f B 2 ) - (D- l)J F \m(A,B)\ 2 -h A (2A 2 + 2A + 2B*B) (2.4) 

The new feature in SU (3)is the occurrence of a term proportional to Jp (B 3 +fi* 3 ) . Such a term is 
permitted by Z(3) symmetry and makes the deconfining transition first order. 
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Figure 2: String tension ratios for thermal gluon potential. 
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Figure 3: Effective Masses for thermal gluon potential. 

3. Potentials 

Particle loop effects contribute to the potential terms in the effective Hamiltonian. The one- 
loop contribution of the massless gauge bosons to V is 

-lf,- A X A {P n ) (3.1) 

B=l U 

with A = 2T 4 /n 2 . This potential is minimized at the elements of Z(N), and yields the familiar 
black-body formula for the free energy density at those points. We have calculated the string 
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Figure 4: String tension ratios for W at B = 1 . 




L 

10 20 30 40 



Figure 5: String tension ratios for W at B = 10. 

tensions for this potential by diagonalizing the quantum-mechanical SU (3) Hamiltonian on a finite 
basis, using as a basis all representations up to eight columns wide in the Young tableaux. In figure 
2 we compare the sorted string tension ratios Or/ Op for Xj K = 1 with Casimir scaling. Although 
there are systematic deviations, they are most apparent for larger representations. We can define 
an effective mass associated with the exponential decay of each representation as a function of 
the scaled time % = (E\ — Eo)t. Figure 3 shows substantial mixing between representations, as 
evidenced by the rapid approach of effective masses to the lowest mass. 

Topological excitations also contribute to the potential. For example, the one-loop calculation 
in SU (2) of caloron contributions to the functional integral in SU (2) at finite temperature indicates 
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Figure 7: Effective Mass Ratios at b = 10 



an instability of the deconfined plasma phase at sufficiently low temperature []1 In certain su- 
persymmetric SU (N) gauge theories on S l x R 3 , confinement can be understoood as arising from 



magnetic monopole effects [11, 12]. The superpotential W, regarded as a function of the Polyakov 
loop eigenvalues, is the afflne Toda potential, and is minimized at Pq. In the case of SU(3), the 
superpotential can be written in the form 



W 



(3.2) 



where b = An/g 2 . The eigenvalues are taken to lie in the first Weyl chamber, where the eigenvalues 
are ordered as Q\ > 02 > 63. We will use W as an example of a potential that is minimized at Pq. 
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Figure 4 shows that the deviation from Casimir scaling increases with the strength of the po- 
tential A as expected. In comparison, figure 5 indicates that at fixed A , Casimir scaling is recovered 
in the limit of large b. As b increases, V is effectively zero except near the boundary between Weyl 
chambers, where two or more Polyakov loop eigenvalues coincide. Because Haar measure leads to 
eigenvalue repulsion, the wave function naturally vanishes when two eigenvalues coincide, and for 
large b, W has little effect on the wave functions. 

The effective masses, shown in figures 6 and 7, show longer plateaus at larger b. This is par- 
ticularly noticeable for smaller representations. This is consistent with the wave functions for large 
b approaching the V = wave functions. Thus increasing b decreases mixing between representa- 
tions of the same 7V-ality. 
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